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Abstract 

We present a generalization of the perturbative construction of the metric operator for non- 
Hermitian Hamiltonians with more than one perturbation parameter. We use this method 
to study the non-Hermitian scattering Hamiltonian: H = p 2 /2m + C-^( x + a ) + C+^( x — a )> 
where £± and a are respectively complex and real parameters and <5(x) is the Dirac delta 
function. For regions in the space of coupling constants £± where H is quasi-Hermitian and 
there are no complex bound states or spectral singularities, we construct a (positive-definite) 
metric operator rj and the corresponding equivalent Hermitian Hamiltonian h. rj turns out to 
be a (perturbatively) bounded operator for the cases that the imaginary part of the coupling 
constants have opposite sign, 3 : (C+) = — 9f(C— )• This in particular contains the 'P7~-symmetric 
case: C+ = C-- We also calculate the energy expectation values for certain Gaussian wave 
packets to study the nonlocal nature of h or equivalently the non-Hermitian nature of H. We 
show that these physical quantities are not directly sensitive to the presence of "P7~-symmetry. 
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1 Introduction 



A Hamiltonian operator H is called "PT-symmetric if it has parity-time reversal symmetry, i.e., 
[H,VT] = 0. Since the publication of the pioneering work of Bender and Boettecher [1] non- 
Hermitian "PT-symmetric Hamiltonians have received much attention. This has led to the discovery 
of a number of interesting theoretical [2j [3j HJ El El [7J as well as experimental [8] implications of 
PT-symmetric Hamiltonians. For extensive reviews see [91 [10] and references therein. 

Perhaps the most prominent feature of a non-Hermitian PT-symmetric Hamiltonian H is that 
its spectrum is symmetric about the real axis of the complex plane. In particular, if H has a discrete 
spectrum, either it is purely real or the nonreal eigenvalues come in complex-conjugate pairs [U [7J. 
It turns out that this is a characteristic property of a wider class of non-Hermitian Hamiltonian 
operators called the pseudo- Hermit ian operators [III [El [13] . A Hamiltonian H is said to be pseudo- 
Hermitian if its adjoint satisfies 

H ] = riHrj' 1 , (1) 

for some Hermitian invertible operator rj. Under the assumption of the diagonalizability of H, one 
can show that its spectrum is real if and only if there exists a positive-definite (metric) operator rj 
satisfying the above equation [TUl [121 US] . In this case H is called a quasi- Hermitian operator [14J . 

The diagonalizability of an operator is equivalent to the lack of exceptional points and spectral 
singularities [15J. Exceptional points are degeneracy points where some of the eigenvectors of the op- 
erator coalesce. This phenomenon have been the subject of many theoretical (HQS] and experimental 
[T7j studies. It may appear for operators acting in finite or infinite dimensional Hilbert spaces. In 
contrast, spectral singularities can only appear for non-Hermitian operators whose spectrum includes 
a real continuous part (See [H] and references therein). Mathematically, they are responsible for a 
break down of eigenfunction expansion [15]. Physically, they correspond to resonances having a real 
energy (zero width) [31 El [19] . 

As we mentioned above, a quasi-Hermitian Hamiltonian is a diagonalizable operator with a com- 
pletely real spectrum. This is not however sufficient reason for using quasi-Hermitian operators as 
observables in a quantum theory. This is because the diagonalizability of an operator and the reality 
of its spectrum do not necessarily imply the reality of the expectation values of the operator. The 
latter condition is in fact equivalent to the Hermiticity of the operator [10] . The advantage of quasi- 
Hermitian operators over other non-Hermitian operators is that they can be made Hermitian by an 
appropriate modification of the inner product on the Hilbert space. This is done using positive- 
definite metric operators rj that satisfy ([1]). The modified inner product is given by (-|-)^ := (-|^-), 
where (-|-) is the inner product that defines the original Hilbert space %. Endowing the vector space 
of state vectors with the inner product (-|-)^, we define a new Hilbert space % p hy S in which H acts 
as a Hermitian operator fTOj, [2Q] . Hereafter we assume that if is a quasi-Hermitian operator and call 
"Hphys the physical Hilbert space. 

In general, rj is not unique. This means that either one must choose rj directly or fix it indi- 
rectly by demanding that a so-called compatible irreducible set of quasi-Hermitian operators will 



act as Hermitian operators in H p h ys , [HI- As explained in [10], the latter approach is very difficult 
to implement in practice, if the only available information is the form of the quasi- Hermitian oper- 
ator H . This is because to select the members of a compatible irreducible set of quasi-Hermitian 
operators containing H, we need to construct the most general metric operator rj fulfilling ([1]). For 
a quasi-Hermitian Schrodinger operator H = — 4^ + v(x) with a typical complex potential, this is 
an extremely difficult open problem. Dealing with this problem is particularly difficult for complex 
scattering potentials such as the one studied in the present article, because the continuous spectrum 
of H is doubly degenerate. In what follows, we will assume that a choice for rj and consequently 
Hphys is made a priori. 

Because rj is positive-definite, it has a unique positive-definite square root p := ^ffj : H — > H. It 
is easy to show that p : "H p h ys — > H is a unitary operator, i.e., 



It establishes the unitary equivalence of the (Hilbert space, Hamiltonian) pairs: ("H p h ys , H) and (%, h) 
where h := pHp^ 1 , [201 El]. The operator h is a Hermitian operator acting in the original Hilbert 
space %. It is called the equivalent Hermitian Hamiltonian associated with the metric operator r), 
[22] . ("Hphysj H) and (H,h) provide equivalent representations of the same quantum system. They 
are respectively called the pseudo- Hermitian and Hermitian representations [TOj [20] . 

In the pseudo-Hermitian representation we work with the physical Hilbert space ("H p h ys , ('W), 
and the quasi-Hermitian observables H, X := p~ l xp,P := p~ l pp, ■ • • , where x,p, •■ ■ are the usual 
Hermitian observables. In the Hermitian representation we work with the usual Hilbert space H 
and the Hermitian observables h := pHp^ 1 , x,p, ■ ■ ■ . A particle that is described by the state vector 
\ip) € H p h y s and the Hamiltonian H can also described by the state vector p\ip) G H and the 
Hamiltonian h. [101120]. 

Working with the Hermitian representation has the advantage of revealing the underlying classical 
system. This is of great importance to derive the physical meaning of the system [23| 121] and establish 
a classical-to-quantum correspondence principle. In order to employ the Hermitian representation, 
we need to compute the equivalent Hermitian Hamiltonian h. This in turn requires the calculation of 
p. A well-known method of constructing h is to use the exponential representation 77 = e~® for the 
metric operator and apply the perturbation scheme developed in [23, 25J. We will begin our analysis 
by extending this method for the cases that the Hamiltonian H involves more than one perturbation 
parameter. We will then apply this method to treat the quantum system defined by the double-delta 
function potential: 



The spectral properties of this and analogous complex point interaction potentials have been consid- 
ered in [26]. See also [271 EH] • A thorough investigation of (jSJ) that addresses the issue of the presence 
of spectral singularities and provided means for locating the regions in the space M of coupling 
constants where the Hamiltonian is quasi-Hermitian is conducted in [T5] . 

In the present paper we offer an explicit construction of an appropriate metric operator for H 
in a three-dimensional subspace of M that includes the "PT-symmetric potentials of the type ([3]) 



(P- IP') = (-\-)v 



(2) 



v(x) = z_5(x + a) + z + 5(x — a), z± 6 C, a £ M + . 



(3) 



with z_ = z*,. Using this metric operator we determine the equivalent Hermitian Hamiltonians 
and compute energy expectation values for some Gaussian wave packets. Our results are valid 
irrespective of the presence of "PT-symmetry. Therefore, they allow for a direct examination of the 
physical consequences of imposing PT-symmetry. 



2 Spectral Representation of Metric Operator for a Scat- 
tering Potential 

Consider a non-Hermitian Hamiltonian H with a complex-valued scattering potential v(x) : R — > C. 
Suppose that v(x) depends linearly on a set z := (z±,Z2, ■ • • ,Zd) of complex coupling constants so 
that iT can be obtained by replacing z with z* in v (x). A typical example is the double-delta 
function potential ([3]). For the cases that H has no bound states (square-integrable eigenfunctions), 
its spectrum is doubly degenerate and its eigenvalue equation takes the form 

H\< k ) = kV a , k )- (4) 

Here a G {1,2} is the degeneracy label and k e R + is the spectral label. In the absence of spectral 
singularities if is a diagonalizable operator, and we can use |Vofc) together with a set of (generalized) 
eigenvectors \4>l k ) of H* to construct a complete biorthonormal system, i.e., 

(€MJ = 5 ab 5(k-q), K q )(€, k \dk = l. (5) 

In this case the following formula gives a (positive-definite) metric operator for the Hamiltonian H 



n 



0=1 



E / \<k)(€, k \dk. (6) 



It is easy to see that rj satisfies (JXJ) We can use this operator to define the positive-definite inner 
product (-|-)^ := (-|^-) and the corresponding Hilbert space H p h ys in which H acts as a Hermitian 
operator. 

Next, we recall that because of the arbitrariness in the choice of the biorthonormal system (in 
particular \<f>* k )), the metric operator (jHJ) is not unique [30]. In what follows we will try to choose 
the biorthonormal system k ), \<f>% k )} in such a way that in the Hermitian limit, where all the 
coupling constants are real, the metric operator tends to unity. 

Since W can be found by replacing z with z* in H, the simplest way to find a set of eigenvectors 
of iff is to replace iTwith z* in \ipl h ). It is however not difficult to see that in general {jV'ofc}) iV'ffc)} 



1 Here we assume that H has no bound sates. If there are N bound sates with real energies, (JB]) is generalized 

il 



t0 V = Eoii J2f=i \<t>tj){<t>ij \ + Eo=i Jq 00 l^lfc)(^f,fcl dft ' where a 3 denotes to the degree of degeneracy of the j-th 



eigenvalue. 



does not satisfy (jSJ). In fact, one can show that 

where X = (Jf t,) is a 2 x 2 matrix that is generally different from the identity matrix. 

One way of constructing a metric operator with appropriate Hermitian limit (namely identity 
operator) is to find a matrix U(z\ k) satisfying 

U\z*-k)K(z-k)U(z-k) = I 2x2 (8) 

and use the biorthonormal system {\4>l k ), \4>lk)} defined by 

:= U ab (z*; k)WQ, \r; jk ) := U ab (z; k)\^ k ) 

This approach relies on the solution of Eq. (jSJ). It is clear form ([7]) that 

K\z*;k) = K{z;k). (10) 

In view of this identity we can rewrite (jHJ) as 

W f (r ; fc) i^i*; fc)s fe)3 ZY(z; fc) = J 2x2 . (11) 

This equation has infinitely many solutions. Perhaps the simplest solution is 

U = K~*. (12) 

In general, has an extremely complicated form. This leads to serious computation difficulties 
in the perturbative expansion of the metric operator. Furthermore, there is no assurance that this 
choice of the biorthonormal system yields a bounded metric operator. 

In Ref. [31] this method is employed to calculate a metric operator for a delta function potential, 
v(x) = z5(x), with a complex coupling constant z having a positive real part. In this case, it yields 
a perturbatively bounded metric operator. We will discuss the possibility of applying this method 
for the complex double-delta function potential (j3J) in Section 4. 



5{k-q)K 



(7) 



(9) 



3 The Perturbative Expansion of h 

Let H : % — > % be a quasi- Hermitian Hamiltonian of the form 

d 

H = HW + HM, H^:=Y,^Hi (13) 

i=l 

where H^°',Hi,-- - , Hd are Hermitian operators and z\, z 2 , ■ ■ • , Zd G C are complex parameters. 
Suppose that \zi\ 1, for all i G {1, 2, ■ • • , d}, so that we can use them as perturbation parameters. 

E 



Then, f lT3|) is a perturbative expansion of H, with and respectively denoting the zeroth 
and the first order terms. 

Consider the perturbative expansion of an arbitrary operator A depending on zi, Z2, ■ ■ ■ , Zd- Let 
rii, n,2, ■ ■ ■ rid be non-negative integers and N := ri\ + ri2 + • • • + rid- Then we call the sum of terms 
proportional to z^z^ 2 • • • z n d d "the iV-th order term" of this expansion and denote it by A^ N \ We 
also use 0(z ) to label the sum of the terms of order greater than or equal to N. 

Because the first order term of the Hamiltonian is generally non- Hermit ian, we write it as 

H V=h£> + h21, (14) 

where H^' and stand for the Hermitian and anti-Hermitian parts of respectively. In view 
of quasi-Hermiticity of H, there is a positive-definite metric operator 77 satisfying W = tjHt]' 1 . This 
relation together with Eqs. ( |T3|) and ( |T4l) imply 

rjHr]- 1 = # (0) + - H®. (15) 

Our aim is to use this equation to construct a perturbative expansion for a metric operator with 
correct Hermitian limit: 

v = i + r] ^ + v {2) + 0{z 3 ), (16) 

and the corresponding equivalent Hermitian Hamiltonian h. 

First, we recall that because rj is a positive-definite operator, there is a Hermitian operator Q 
satisfying [25] 

rj = e' Q . (17) 

Next, we use (fT5j) . ffTT|) . the Baker- Campbell-Hausdorff identity, and the perturbative expansion of 
Q, namely 

Q = Q^ + Q^ + 0(z 3 ), (18) 

to obtain EDI: 



Ft = H+[H,Q] + ^[[H,Q},Q} + 0(z 3 ) (19) 



H<n+HV) + [# (0) ,Q (1) ] + [# (0) ,Q (2) ] + [# (1) ,Q (1) ] + \[[HV>,QV)],QM]+0{2?). 



Comparing the right hand side of (fl5|) with the first order term of the last equation, we have 

hV-hV=hV + hV+[h®qV], (20) 

or equivalent ly 

[H&>,QV] = -2h£L. (21) 



It is also clear from f lT5|) that the second order term of vanishes. This implies that the second 
order term in ffT9l) must satisfy 

[J5T« Q«] + [# (0) ,g (2) ] + l[[H(°\Q(%QU] = 0. (22) 

Similarly, using 

p = e~ Q ' 2 (23) 
we find the perturbative expansion of the equivalent Hermitian Hamiltonian h: 

h = pHp- 1 = H +[H^]+ 1 -[[H,%^] + 0(z 3 ) (24) 

= h^ + hu + 1 i[H^,Q^\ + [# (1) ,g (1) ] + [# (0) ,q (2) ] + ^[[# (0) ,g (1) ],g (1) ]} + o(z 3 ). 

In light of (l2~Tj) and (!22j) . we can simplify this expression as follows. 

= #C°> + F« - J [[#^ Q^ 1 (25) 
o 4 

This is a straightforward generalization of the results of Refs. [23] to the cases involving more than 
one complex perturbation parameter. See also |32j . 

Next, we use the identity Q« = -77M to express the integral kernel h(x,y) := (x\h\y) of the 
equivalent Hermitian Hamiltonian. This yields 

h(x,y) = ^°)(x,y) + 4 1) (x,y)-i(x|[i?W,r ? W]|y) + 0(^) (26) 

= hQ\x,v) + h£\x,v)-\I (^i(x,^ (1 Hy^^)-^ (1 H^yO^^(y , ^))^ / + o(^ 3 )• 



Having obtained the equivalent Hermitian Hamiltonian, we can examine the physical content of 
the model using its Hermitian representation. Alternatively, we can purse the study of this system 
in its pseudo-Hermitian representation. This requires the construction of the pseudo-Hermitian 
observables O = p^ x op, where o = x,p, ■ ■ ■ are the usual Hermitian observables. Following a similar 
approach to the one leading to (12"4"I) . we find 

O = p- 1 op = o-[o,f] + i[[o,f],f] + 0(z 3 ) 

= o- l - |[o,g«] + [o,qW] - J[[ ,gW],g«]} + o(z 3 ). (27) 

4 The Double-Delta Function Potential 

4.1 Eigenfunctions and the A'- matrix 

In this subsection we summarize some of the properties of the double-delta function potential that 
are reported in Ref. [15]. 



Consider the time-independent Schrodinger equation: 

h 2 d 2 



Rip 



(+5(x — a) + C-<K X + a ) 



?p = Eip. (28) 



2m dx 2 

Let I be an arbitrary length scale. Defining the dimensionless quantities 

2m£(± x a 2m£ 2 E , . 

*±--=-f^, v-= T «:= 7 , E --=-V~ ( 29 ) 

we express the corresponding dimensionless Hamiltonian as 

2m£ 2 R d 2 

H := — — = — — + z+5(x - a) + zS(x + a). 30 
nr dx 1 

We can easily solve the eigenvalue problem for H and obtain the following expression for the eigen- 
vectors ipl k of the Hamiltonian H. 



?Pl k (x) = -^{e ikx - ^ [e" ifc(x+2a) - e ikx ] 6{-x - a) - ^± [ e lkx - e -**0»-*O] 8(x - a)), (31) 



V>f» = (32) 
where ^(x) := [sign (a;) + l]/2 is the step function. Clearly, 

^ k (x) = ~^={e ikx ~ ^ [ e - ik{x+2a) - e ikx ] 9(-x - a) - ^ [e lfca; - e - 4fc ^- 2a )] 0(x - a)}, (33) 



VSb(*) = (34) 
and the entries of the matrix K of Eq. ([7]) takes the form 



2 2 

K U (A;) = K 22 (A;) = 1 + ^_!±, (35) 

K 12 {k) = {Ak 2 )- 1 [iz„{2k - iz_)e 2iak - iz+(2k + iz + )e- 2iak ] , (36) 
K 21 (k) = K X2 {-k). (37) 

As we can infer from these equations, the matrix K~ l l 2 has a rather complicated expression. This 
makes a direct application of the method of Section 2 extremely difficult. In what follows we will 
try to pursue a different approach for constructing a metric operator with a correct Hermitian limit. 
We will also demand that, at least to the first order of perturbation, the metric operator be densely- 
defined and bounded [TO] . 

Similarly to the case of the single delta function potential, v(x) = z8(x), the fact that we have 
exact and closed-form expressions for the eigenf unctions of the Hamiltonian ( 1301) does mean that we 
can perform an exact calculation of a metric operator. Note that the latter requires choosing an 
appropriate set of eigenfunctions (f> z ak and evaluating the integral in (JBJ) exactly. Lack of a systematic 
method of selecting <p z ak (alternatively the matrix- valued function U appearing in (jSJ)) is the main 
reason why we conduct a perturbative analysis of the problem. 



4.2 Perturbative Calculation of the Metric Operator 



For the cases that 3l(z±) > and |e±| := 



S(*±) 



<C 1, the Hamiltonian if is a quasi-Hermitian 



operator [15]. Therefore e± := ^r^j may be employed as perturbation parameters. As shown in 
Appendix A, this choice leads to a metric operator that does not tend to the identity operator in 
the Hermitian limit (e± — » 0). An alternative choice for perturbation parameters is the coupling 
constants (z+,Z-) =: z themselves. In the remainder of this section we construct a metric operator 
r] using these perturbation parameters. 

According to (|35|) - (137|) . the zeroth order term of the matrix K obtained for the eigenvectors (}3~Tj) 
- ( 1341) is the identity matrix. This in turn implies that the zeroth order term for the corresponding 
metric operator is the identity operator. Yet this metric operator is plagued with the problem of 
unboundedness and the lack of a dense domain. 

In order to construct a densely-defined and bounded metric operator, we use the following ansatz 
for the eigenvectors of W . 

2 2 
6=1 6=1 

where u±^(k)'s are free weight functions. In Appendix B we describe a procedure for selecting a 
proper set of weight functions. We could do this successfully only for the special case that the 
coupling constants differ by a sign: z + = —z_ =: z. In this case, we find (see Appendix B) 

0i» = (l ~ Q\ + ^T^ "^ + ° {z2) > (39) 

4A X ) = (40) 
T](x,y) = 6(x-y) + i 1 (1 \x,y) + 0(z 2 ), (41) 

V (1) (x,y) = i -^^sign{x - y)[6{x + + y + ) - B{x~ +y~)} (42) 

where 

x := x ± a, y ± := y ± a. (43) 
The metric operator ( |4T|) has the following desirable properties. 

1. It tends to the identity operator in the Hermitian limit. 

2. It is a densely-defined bounded operator. 

3. It satisfies the differential equation for the (pseudo-) metric operators associated with pseudo- 
Hermitian Schrodinger operators |33j. 

4. For the "PT-symmetric case corresponding to a purely imaginary z, it reduces to the metric 
operator obtained in [B]. 



r 



r 



Figure 1: The figure on left shows the curves in the az-plane along which H has 
spectral singularities. Here r := a^ft(z) and s := aQ(z). The vertical dashed lines 
correspond to the PT-symmetric case with purely imaginary z. The diagonal lines 
are the lines s = ±r. The figure on the right shows the number of bound states 
(generally complex eigenvalues with square- integrable eigenfunctions). As the color 
changes from the darkest to lightest the number of bound sates will changes from 
zero to four. 

We would like to emphasize that the above construction is valid only for the cases that the Hamil- 
tonian is quasi-Hermitian. Otherwise the metric operator will not satisfy the pseudo-Hermiticity 
relation ([1]). Therefore, it is of outmost importance to determine the range of valued of z for which 
the Hamiltonian is quasi-Hermitian. These are the regions where H has no spectral singularities 
or complex eigenvalues. Figure Q] shows the regions in the complex az-plane where the Hamilto- 
nian has spectral singularities and bound states. This figure is obtained using the contour integral 
method described in Ref. [15]. For small values of \z\ with > H has a bound state 

(an eigenvalue with a square- integrable eigenfunction). This corresponds to a real eigenvalue if and 
only if \Q(z) \ = 0, [15]. In the darkest region shown in the right-hand figure the Hamiltonian is free 
of spectral singularities and bound states. This is a region where it is quasi-Hermitian, and ( 14TT) 
provides a reliable metric operator. 



i n 



4.3 Equivalent Hermitian Hamiltonian 

Inserting (j4"2|) in ( )26|) and doing the necessary calculations, we obtain the following expression for 
the equivalent Hermitian operator h defined by the metric operator (j4Tj) . 

h(x,y) = 5(x-y)(^-^ + ft(z)[5(x-)-5(x + )}^ (44) 
^ f 5(x+)[0(y+) - 6{y - 3a)] + 5(x~)[9(y + 3a) - 0(y-)] 



+5(y + )[6(x + ) - 9{x - 3a)] + 5{y-)[0{x + 3a) - 0(ar)] J + 0{z 3 ). 

If we multiply h(x,y) by and use (T2§j) . we obtain the dimensionful Hermitian Hamiltonian 
operator: 

h(x,y) = 6( X -y)(-^^ + M(Olt(x-a)-6( X +a))\ (45) 
+ m|C)^| 5(x + + a) _ 6{j _ 3a)] + 5(x _ Q)[0(y + 3a) _ 6{j _ a)] 

+5(y + a)[6(x + a) - 0(x - 3a)] + 5(y - a)[6(x + 3a) -0(x-ot)]} + C(C 3 ), 

where C = C+ = ^j- The second order (nonlocal) part of the Hamiltonian ( 145]) is equivalent to the 
anti- Hermitian part of the Hamiltonian (130|) . In the following we study the effect of this nonlocal part 
on the energy expectation value, := (^Ih^), for a particle described by a normalized Gaussian 
position wave function ip G L 2 (IR). 

The action of h on an arbitrary element ij; of the Hilbert space T-L is given by 

h 2 

(x|h|z/>) = f fx) + R(C) [i(x - a)^(a) - i(x + o)^(-a)] (46) 

2m 



+ Ap$ 1 ^ x + «) - #( x - 3a)]^(-a) + [0(x + 3a) - 0(x - a)]^(a) 



m[9(C)]_ 

^(y)dy + 5(x- a) / 

-a J— 3 



3a 

^(y)dy I + 0(C 



The first line of this Equation coincides with the action of the Hermitian part of the original quasi- 
Hermitian Hamiltonian, namely H^. = H^ ) + H^. 
In view of (l4"6l. 



fc2 poo 

= —J W(x)\ 2 dx + K(()[ma)\ 2 -m-a)\ 2 ] (47) 



-oo 

m[S(C)] 2 ' 



2h 2 



3u 



F "" L ^ /J 3ft(V*(-Qf) / ^(x)dx + ^*(a) / ^(x)dx ) + C(C 3 ) 



where ^ is a normalized wave function. A typical example is a Gaussian wave packet, 



with mean position (x)^ := (ipi\x\ipi) = and mean momentum (p)^ := (^iIpI^i) = hk. Substi- 
tuting gH]) in (dZD, we find 



% = ^ (2A: 2 - o~*) + U(a, a, k) + 0( C 3 ), 



where 



?7(a, a, k) := exp 



a 2 + k 2 a 4 
2a 2 



-ika 



. ika 2 + 3a\ / ika 2 — a 

erf = — erf 



y/2. 



a 



V2a 



(49) 



(50) 



describes the effect of the nonlocal part of h (equivalently the non-Hermitian part of H), and erf (x) := 
^= J* e~ t2 dt is the error function. Figure [2] shows the plots of U(a, a, k) and U(a, a, 0) for a = 1. It 




u 

1.20r 



LO 1,2 1,4 1,6 1.8 2.0 



Figure 2: Plots of U(l, a, k) (on the left) and £7(1, a, 0) (on the right). 



turns out that the non-Hermiticity effect attains its maximum around (a, k) = (1.5a, 0) and decays 
rapidly for the mean momentums (p)^ = hk outside the range (— ^, ^). 

Next, we compute the energy expectation value for a stationary Gaussian wave packet with mean 
position (x)^, 2 = x : 

\2 N 



^ 2 (x) 



1 



, /TO" 



2)1/4 



exp 



ix - Xpj 
2a 2 



In view of (1471). we have 



E„ 



h 2 



>P2 



+ 



V(a, a, x ) + -L-M-W(a, a, x ) + 0(C 3 ) 



V(a,a,x ) : 
W(a,a, x ) := exp 



4mcx 2 a^ v ' ' 2y/2tP 
exp (— er~ 2 (x — a) 2 ) — exp (— (x~ 2 (x + a) 2 ) 



(a + x ; 
2cr 2 



erf 



+ e" 



a + x 
y/2> 



a 



erf 



+ erf 

a — x 



^2, 



(7 



3a — x 
+ erf 



3a + x 
y/2, 



a 



Here W(a, a, xq) reflects the effect of the nonlocal part of h (non-Hermitian part of H). 



(51) 



(52) 
(53) 
(54) 



1 o 



Figure [3] shows the plots of W(l, a, Xq) for o e [0.03,0.6] and a £ [0.2,10]. The non-Hermitian 
effect attains its maximum at (xq ~ 0, a ~ 1.5a). For small values of a, it persists for wave packets 
with mean position belonging to open intervals (±a — 1.5a, ±a + 1.5a). Outside these intervals it 
decays rapidly. 




Figure 3: Plots of W(a,a,x ) for a e [0.03,0.6] (on the left) and a e [0.2,10] (on 
the right) in units where a = 1. 



4.4 Pseudo-Hermitian Position and Momentum Operators 

To calculate the dimensionless pseudo-Hermitian position and momentum operators corresponding 
to the metric operator (f4TT)-( f42|) . we substitute x and p for o in ( 1271) and use = — r/ 1 ) and the 
identities: 

(x\[x,A]\y) = (x-y)(x\A\y), (x\\p,A]\y) = -i{d x - d y )(x\A\y), 
where A is a linear operator. This yields 

(x\X\y) = ( x \x\y) + l ^\x-y\[9(x + + y + )-9(x-+y-)]+0(z 2 ), (55) 
(x\P\y) = (x\p\y} + %(z)5(x-y)[6(x + )-6(x-)]+0(z 2 ). (56) 

Next we obtain the dimensionful pseudo-Hermitian position (X := £X) and momentum (P := |P) 
operators: 

(x|X|y) = ( x | x | y ) + !!^| x -y|[^(x + y + 2a)-^(x + Z/ -2a)] + 0(C 2 ), (57) 
(x|P|y) = ( x |p|y) + _^5( x - y )[0( x + a)-0( x -«)] + O(C 2 ). (58) 

4.5 Calculating Metric for More General Cases 

In Section 4.2, we constructed a metric operator with the desired Hermitian limit for the cases that 
the coupling constants (± differed by a sign. Our construction was based on the spectral method 



that yielded r\ in terms of its spectral decomposition. An alternative method of constructing a metric 

2m dx 2 



operator for a Schrodinger operator, — ^- + v(x), is the one based on the universal differential 



equation [33] : 

2m \ 
-dl + d 2 y + w K(x) - „(y)]J V (x,y) = 0. (59) 

In this section we use this method to extend the results of the preceding sections to a more general 
class of double-delta function potentials. 

Consider a quasi- Hermit ian Hamiltonian operator, 

H^HW+HW =-i + t^(x), (60) 

with a purely imaginary potential v% = i^s(vi), and a corresponding metric operator r]i satisfying 
( 1591) with v = v\. Let ( be a complex perturbation parameter such that t>i is proportional to 3(C). 
This suggests the following perturbative expansion for 77. 

Vl = l + V ^ + 0(C 2 )- (61) 

Next, suppose that the potential v± is supplemented with a real part v% that is proportional to 
^ft(C)- Then it is easy to see that up to the first order of perturbation, 771 satisfies fl59|) for the potential 
i>i(x) + f 2 (x), i.e., r]i is a metric operator also for the Hamiltonian 

H 2 = Hi + h£> = + Wl (x) + t; 2 (x). (62) 

Furthermore, in light of f l25|) . the equivalent Hermitian Hamiltonian for Hi and H 2 are respectively 
given by 



hi = H(°)-i[Hii,^] + (9(C 3 



h 2 = H(°) + H«-i[Hii,r / «] + 0(C 3 )=h 1 + H« + 0(C 3 ). (63) 

Now, we confine our attention to the double-delta function potential. In Section 4.2, we con- 
structed an appropriate metric operator, namely for a double-delta function potential whose 
couplings differed by a sign. In view of the argument given in the preceding paragraph, it is also 
a valid metric operator for the more general case that the real part of the coupling constants are 
arbitrary (but small) and their imaginary part differ by a sign: 

H = -|l + C+5(x-a) + C-5(x + a), 3(C+) = -3(C-)- (64) 
2m 



Another class of potentials that admit ( f4"Ti) as an appropriate metric operator is [23 



v(x) = -£5(x)+i\{5(x-a)-5(x + a)}, £ £ IR + , A G R. (65) 



1 .1 



Next, we use the metric operator (j4ip to compute the equivalent Hermitian Hamiltonian h for 
the Hamiltonian f )64l) . Using ( |63i) and performing the necessary calculations, we find 



As we expected the nonlocal part of h is identical with the one obtained for the special case where 
coupling constants differ by a sign. 

The general PT-symmetric double-delta function potential (C+ = C* =: C) corresponds to a 
special case of the Hamiltonian Our analysis shows that, within the framework of perturbation 
theory, the physical effects associated with the non- Hermit icity of the Hamiltonian H, that are 
contained in the nonlocal part of the equivalent Hermitian Hamiltonian h, are not sensitive to the 
presence of "PT-symmetry. This is because by perturbing the Hermitian part of the Hamiltonian we 
may destroy its PT-symmetry while preserving the same non-Hermitian (nonlocal) effects on the 
physically observables quantities like energy expectation values. 

Finally, we wish to stress that the above calculations of the metric operator and the equivalent 
Hermitian Hamiltonian are reliable only within the regions in the space of coupling constants that 
the Hamiltonian does not posses spectral singularities or bound states with complex eigenvalues. In 
Figure H] we give the location of the spectral singularities and the number of bound sates for the 
general "PT-symmetric casein In the region painted by the darkest color in the r-s-plane (where 
t '■= ^|r^9ft(C)> s := npr^tO) the Hamiltonian has no spectral singularities or bound states. Hence 
in this region it is quasi-Hermitian, and ( HTT) gives a reliable metric operator provided that we stay 
within the part of this region that is close to the origin. 

5 Concluding Remarks 

In this article we have employed the pseudo-Hermitian formulation of quantum mechanics to study 
a quantum system defined by a Hamiltonian with two complex point interactions, H = p 2 /2m + 
(-5(x + a) + C+<K X — a )- This requires the construction of an appropriate metric operator that 
reveals the structure of the physical Hilbert space and also the observables of the theory. It further 
allows for the construction of an equivalent Hermitian Hamiltonian operator. 

The main difficulty one encounters in trying to construct a metric operator for H is that depending 
on the choice of the eigenfunctions of W, one obtains a "metric operator" that may be ill-defined 
or unbounded. In this article we could successfully construct a densely-defined and bounded metric 
operator to the first order of perturbation for the special cases that the coupling constants £± differed 
by a sign. We use this metric operator to compute the corresponding equivalent Hermitian operator. 

2 Again we have used the contour integral method described in Ref. [TS] to obtain the location of the spectral 
methods and the number of bound states. 




(66) 
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Figure 4: The curves along which the PT-symmetric Hamiltonian ( 164)) with ( + = 
C* =: ( has spectral singularities (the figure on the left) and the number of bound 
states with (the figure on the right). As color changes from the darkest to the 
lightest, the number of bound sates takes the values 0,1,2,4, respectively. The 
vertical dashed lines correspond to the purely imaginary coupling (, and the diagonal 
(dotted-dashed) lines are the lines r = ±s. There is a bound state within the circle 
(r + l/2) 2 + s 2 = l/4. 

This in turn allowed us to compute energy expectation values for a class of Gaussian wave packets. 
We then investigated the physical consequences of the non-Hermiticity of the Hamiltonian H by 
examining the contribution of the anti-Hermitian part of the Hamiltonian (equivalently the nonlocal 
part of the equivalent Hermitian Hamiltonian) to the energy expectation values. 

In view of the fact that the integral kernel of the metric operator is a solution of a certain 
differential equation, we could generalize our results to the cases that only the imaginary part of the 
coupling constants differed by a sign. This allowed for the application of our results for a general class 
of double-delta function potentials that included all PT-symmetric double-delta function potentials 
as a subclass. Our investigation of the physical effects of non-Hermiticity shows that (to the first 
nontrivial order of perturbation theory) these effects are not directly sensitive to the presence of PT- 
symmetry. This is because we can easily perturb the real part of the potential in such a way that 
the effects of non-Hermiticity is left unaltered while PT-symmetry is destroyed. Note however that 
if such a perturbation does not violate the quasi- Hermit icity of the Hamiltonian, the Hamiltonian 
will necessarily possess a symmetry that similarly to the PT-symmetry is generated by an invertible 
antilinear operator (T3J [30] . This symmetry cannot however be interpreted as the space-time reflection 
symmetry. 
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Appendix A: Another Choice for the Wight Functions 



if 



3(*±) 
RCb±) 



<C 1, and 3?(z±) > 0, the Hamiltonian (130]) is a quasi- Hermit ian operator [15], and we can 



use e± := as perturbation parameters for a perturbative calculation of a metric operator. This 
requires selecting an appropriate set of \4> z ak ) that would lead to a densely-defined bounded metric 
operator. A natural candidate is the following direct generalization of the expression obtained for a 
delta function potential with a complex coupling constant [3Tj . 

_i 

\<Pl k ) = u(z+,Z- ] kM; k ), u {z +1 z^k):=[l+ Z -^ \ (67) 



where 7 G M. + is arbitrary and \iftf k ) are given in ( )33|) - ([3 
In view of ®, (JMJ) and (JB7D, 

/■oo 

7 7 (x, 3 /)= / < fc (x)« fc (|/))*IU(z + ,z_,A;)^, z_, fc) := z., k)\ 2 . (68) 



Introducing k := -^== '■= pk, E\ := e + + e_, and e 2 := e+e_, we can expand W = W(z + , z_, k) 



W (r + (l + ie + ), r_(l + ze_), ^ in powers of e: i 



1 



Employing this relation in ( 1681) and performing the necessary calculations, we find 

rj(x,y) = r]o,o(x,y)+r} +i+ (x,y) + r}-„(x,y) 

+T,+ (x,y)+ri-(x,y)+Ti- + (x,y) + (x yf + C(e 4 ), (70) 



Here stands for the sum of terms proportional to t\^_ with a + b = N . 



1 17 



where 



x-y 
P 



V±,±{x>y) = l r 2 ± {l + e 2 ± )6{±x*)6{±y*) £. 



ja ( x-y \ _ fi ( x T + y T 



p 



v±(x,y) 



r±(l + ie±] 



V „ i+ (x, y) = ±^9(y-)9(-x + )(l + e 2 + - e_)) x 



8(±y T ) S. 



i (x-y 
P 



i f x T + y T 



P 



(71) 
(72) 
(73) 



£. 



^2 ( x + + y + ^ + f2 ( x + y~ 



P 



?2 i x-y 



P 



T 



2l x-y + 4a 



P 



£:= (l,e 2 ,--^) = (l,e+e_, 



X n (2^,1,2^1,2,2^,3), n — 0,1,2, 



1 f 00 fcl—n^ika _^ 
In.m(Oi) := — / ~ T7T7— 0?fc, X m 1= (Xo )7Tt , Xl )m , X2, m ) , 771 > 1, 



Xi = (a 



2tt (1 + p)' 

e -M z ' e -|a| 



(74) 

(75) 
(76) 



-|a| 



-sign(aj 



X, 



-|«| 



X, 



16 



2 ' 2 ° v 7 ' 2 / ' 4 
1 + |a| — a 2 , ia(l + |a|) , 3(l + |a|) + a 2 ), 



(1 — \a\ , za , 1 + |a|) , (77) 

(78) 



x ± and y ± are given by (1431) . and (x -H- y) stands for the sum of the previous terms with x and y 
interchanged. 

It is not difficult to see that £ = (1, 0, 0) + 0(e 2 ) and £ ■ X n = X nil + 0(e 2 ). Inserting the values 
of X n i given by (177)) - (1751) in (17T1) - (1741) and rearranging the terms, we have 

r)(x,y) = \5{x-y) 



1 2 

+' '^l^- + P[rl9(x-)9(y-) + r 2 _9(-x + )9(-y + )] - ^-0(-. r + )tf(//-)(:i -[ ,y_ - ; f _ ) 



+e p 



Up 8' 

+ ^ 1 6'(2/~)sign(a: - y)(l + ie+) - —9{-y + )sign(x - y)(l + ie_)| 
■ ^000%-) + ^0(-x + )0( 2/ -)(l + *e + - *e_) 



% )sign(x +y )(1 +«e+)| 



+e" 



■0(-y + )sign(x + + y + )(l + ie_) J 



-e — I ^0(-x + )0(tT )(1 + ie+ - ie_)j 



+ (x^y)* + e>(e 2 ). 



(79) 



As suggested by this relation, 77 is actually densely-defined and (perturbatively) bounded, but in the 
Hermitian limit it does not tend to the identity operator. 



Appendix B: A Proper Choice for the Wight Functions 

One can check that unless one fixes u±$(k) in a very special way the metric operator corresponding 
to ( 1381) is an unbounded operator. In order to obtain such a special choice we consider the ansatz: 

<%k{*) = ^) + ^) + 0(z 2 ), (80) 

4A X ) = <-*0«0> ( 81 ) 

where ip z a \{x) are given by Eqs. fl33l) - fl34|) . 



A=± j=l 

m 



m 

E E («;,^(*) + ^2,a^^)) + <^ 2 ) 



- E E T Ki^" 1 + ^AA^"* ) + ^), (82) 

A=± j=l 

m is a positive integer, and 6j and Vj t b,\ are respectively real and complex free coefficients. By defining 
fi(b) := 2b — 3, we can rewrite 0° fc (x) in a more compact form, namely 

= 4= E E E l -^e-^^ + 0{z% (83) 

V ZlX A=± b=l j=l 

For the cases where flHTj) holds, we can use (jSJ) to express rj(x,y) = (x\rj\y) as 

/DO 
0jik(*)«*(y))*d*- (84) 
-oo 

Using 

i ;>oo ifca i roc 

— —dk = sign{a), — e ika dk = 5{a), (85) 
and ([HOD " ©, we then find 

V {x, y) = S(x -y) + V V ( x , y ) + 0(z 2 ), (86) 



where 



r] (1 \x,y) := r) { + '(x, y) + rjL '(x, y) + r]^ '(x, y) + (x y)* , (87) 
V+\ x ,y) := ^ {sign(x + y - 2a) - siga{x - y)}6(y - a), (88) 

V-\ x ,v) := -^{sign(x + y + 2a) - sign(x - y)} 6(-y - a), (89) 

2 m 2 m 

r£\x,y) := EEEto^EEE^^ + ^H")- (90) 

A=± b=l j=l A=± 6=1 i=l 

Note that ii^\x,y) depends on our choice of the free coefficients Vj^x and bj in (|83|) . Our aim is 
to make a choice that would lead to a an appropriate metric operator. 



i n 



First, we use the identity 

sign(w + t>)[sign(u) + sign(u)] = 1 + sign(w)sign(t>), (91) 

to obtain 

V+\x,v) = y[-2sign(x - y) - sign(x~ + y~) + 1] + ^sign(x - y)9(-x~ - y~) 

= y[sign(oT + y~) + 1] - ^sign(x - y)9{x~ + y~), (92) 

T)-\x,y) = — [2sign(x — y) — sign(x + + y + ) + 1] — — sign(x — y)9(x + + y + ), (93) 
8 4 



where x and y are given by Eq. (H3|) . 
Next, we rewrite 4>\ k (x) in the form 



where 



Z + -ik(x-2a) -ik(x+2a) , Z - ikx 

4 4 + 2 



(94) 



2 m' • * * 

izXv, 



2vr t , . i • 1 ^ 

A=± b=l 3 = 1 

and m' is a positive integer. Using the above equations, the first order term i]^\x,y) of rj(x,y) can 
be rewritten as 

r/ w (x,y) = rf^(x,y) + v'-\x,y) + rf™(x,y) + (x y)* , (95) 
vi\x,y) = ^[l-2sign(x-y)e(x*+y*)\, (96) 

2 m' 2 m' 

vf\x,y) = E E E & (^. ?/) = E E E ^^(^ + + y ( 97 ) 

A=± 6=1 3=1 A=± 6=1 3=1 

= rlp{x,y) + ysign(aT + y~) - ysign(x + + y + ) + ^sign(x - y). 

Setting (f>'\ k (x) = or equally r]'^'(x,y) = yields 

V ,(1 \x,y) = rf^(x,y)+rf ( l\x,y) + (x o y)*, (98) 

which is generally not a bounded operator. But if we confine our attention to the special case where 
z + = —Z- =: z, we find 

V ' {1 \x, y) = ^sign(a; - y)[6(x + + y + ) - B{x~ + y')], (99) 
which is a bounded operator [6]. 
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